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We present the ETMC results for the bag parameters describing the neutral kaon mixing in the 
Standard Model and beyond and preliminary results for the bag parameters controlling the short 
V~j distance contributions in the D° — Z5° oscillations. We also present preliminary results for the Bg c /, 

Bbs, Bbs/BbcI and E, parameter controlling B Q d ,—B d j s oscillations in the Standard Model employ- 
ing the so-called ratio method. Using Nf = 2 maximally twisted sea quarks and Osterwalder- 
Seiler valence quarks we achieve both 6 (a) -improvement and continuum like renormalization 
pattern. Simulations are performed at three-values of the lattice spacing and several values of 
quark masses in the light, strange, charm region and above charm up to ~ 2.5m c . Our results are 
extrapolated to the continuum limit and extrapolated/interpolated to the physical quark masses. 
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1. Introduction 



Flavour Changing Neutral Currents (FCNC) and CP violation may provide relevant informa- 
tion on the impact of beyond the Standard Model (BSM). The lattice computation of the relevant 
matrix elements appearing in K° - K and B° - B mixing in combination with the experimen- 
tal value of Ek , AM Bc j an d AM^ offers the opportunity to constrain the BSM model parameters. 
D°-D u , at variance with K° - IT and 5° - B mixing, involves valence up-type quarks so it is 
sensitive to a different sector of New Physics (NP). Although long distance effects dominate over 
short distance effects in the D system, it is still possible to put significant constrains on the NP 
parameter space. The most general BSM AF = 2 effective Hamiltonian of dimension-six operators 
contributing to neutral meson mixing is 

H^ 2 = tQ(m (i.i) 

where the operators <2, involving light (/) and non-light (h) quarks read 

Qi = [h a r^-rs)i a ] [h%(i-r 5 )i b ] 

Q 2 =[h a (l-y 5 )l a ][h b (l-Y 5 )l b ] Q3=[h a (l-r 5 )l b }[h b (l-7 5 )l a ] (1.2) 
q 4 = [£»(!- )$)/«] [h b (\ + y 5 )l b ] Q 5 = [h a (l-75)l b ] [h b (l + Ys)l a ] 



2. Lattice setup 



Our lattice computations has been performed at three values of the lattice spacing using the 
Nf = 2 dynamical quark configurations produced by ETMC. In the gauge sector, the tree-level 
Symanzik improved action has been used while the dynamical sea quarks have been regularized 
employing a twisted mass doublet at maximal twist [1] which provides automatic G(a) improve- 
ment [2]. The fermionic action for the light doublet y in the sea reads in the so-called physical 
basis 



n.Mtm 




+v; 



li 



(2.1) 



In addition, both €?(a) improvement and continuum-like renormalization pattern for the four- 
fermion operators are achieved by introducing an Osterwalder-Seiler [3] valence quark action al- 
lowing for a replica of the heavy (h,h') and the light (1,1') quarks [4]. The valence quark action 
reads 



r.OS 

°val 



E E <?/ 

x f=\.V,h.h' 



1 



:rn(Vn + Vl)-ir 5 r f 



Mr 



E V M V M 



+ Hf\q f (x) 



(2.2) 



where the Wilson parameters are conveniently chosen such that r/, = r\ = ry = —ry. 

In table 1 we give the details of the simulation and the values of the sea and the valence quark 

masses at each value of the gauge coupling. The smallest sea quark mass corresponds to a pion 

of about 280 MeV for the case of /3 = 3.90. We simulate three heavy valence quark masses jU<y 

around the physical strange one to allow for a smooth interpolation, further three n> c - around the 

physical charm mass followed by a sequence of heavier masses in the range (m c ,2.5m c ). For the 
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a-\L 3 


XT) 








a\i h 


3.80 


24 3 x 


48 


0.0080 0.0110 


0.0165 0.0200 0.0250 


0.1982 0.2331 0.2742 


0.3225 0.3793 .4461 


(a~().lfm) 












0.5246 0.6170 


3.90 


24 3 x 


48 


0.0040 0.0064 


0.0150 0.0220 0.0270 


0.1828 0.2150 0.2529 


0.2974 0.3498 0.4114 








0.0085 0.0100 






0.4839 0.5691 


{a ~ O.085fm) 


32 3 x 


64 


0.0030 0.0040 








4.05 


32 3 x 


64 


0.0030 0.0060 


0.0120 0.0150 0.0180 


0.1572 0.1849 0.2175 


0.2558 0.3008 0.3538 


{a ~ 0.065fm) 






0.0080 






0.4162 0.4895 



Table 1 : Simulation details 

inversions in the valence sector we used the stochastic method with propagator sources located at 
random timeslices in order to increase the statistical information [5, 6]. Gaussian smeared quark 
fields [7] are used for masses above the physical charm one 1 in order to improve the determination 
of the ground state contribution with respect to the case of simple local interpolating fields. The 
value of the smearing parameters are kc = 4 and No = 30. In addition, we apply APE-smearing to 
the gauge links in the interpolating fields [8] with the parameters GCape = 0.5 and Nape = 20. 

3. Bag parameters 

The bag parameters associated to the operators in Eq.(1.2) are denned as 



{P v \Qi{»)\n=C l B l (p)m P ft = C l B K (ji)m P ft 



(P u \Q i ( y l)\P°)=C i B i (v) 



mpfp 



(3.1) 



where C, = 8/3, —5/3, 1 /3,2,2/3, i = 1, ..,5. The renormalisation constants (RCs) of the relevant 
four- and two-fermion operators has been computed non-perturbatively in the RT-MOM scheme 
[9]. Due to the OS-tm mixed action setup, the renormalized values of the bag parameters are given 
by the formulae [4] [10] [11] 

Zn 



Si, 



B, 



—^-Bj 
ZpZ s 



i,J = 2 ) ..,5 



(3.2) 



4. K°-T 



For large time separation yo <^ xq <C yo + T sep the plateau of the following ratio estimates 

C l(*0) c r Dl/ s Q(X ) 



E\Bi\(xo) 



E[Bi}(x ) 



Cap(jco)C^p(*o)' Cpp(x )C' P p(xo) 
provides the estimate of bare B,-. The involved correlators are 

Q(x ) = l d (^ +Txp Q i (x,x Q )^) 

Cpp(xo) =Z(P 12 (x,x )^), C AP (x ) = £<A 12 (^o)^i> 
C'pp(xo) = i(^ +T P 34 frxo)), C' AP (x ) = £(^ + r A 34 (f,x )) 



(4.1) 



(4.2) 



see the last two columns in table 1 
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MS(2 GeV) 


Si 


B 2 


S 3 


B 4 


B 5 


0.52(02) 


0.54(03) 


0.94(08) 


0.82(05) 


0.63(07) 



Table 2: Continuum limit results for fi; parameters of the K° — K system renormalized in the MS 
scheme of [12] at 2GeV using Ml -type RCs defined in [11] 



(a) (b) 
Figure 1: (a) Comparison of the plateau for the estimator of the B\ bag parameter at j3 = 3.80 on a 
24 3 x 48 lattice and (a/i/,a/X/j) = (0.0080,0.2331) with local and smeared sources (b) Plateaux for 
E[Bi] (i=l,...,5) plotted vs t/T sep for j3 = 4.05 , (a^ h a^ h ) = (0.0030,0.1849) on a 32 3 x 64 lattice. 
The shaded region delimits the plateau interval. 

with 2 &n = Y,q2(y,yo)75qi(y,yo), ^ = ^4(9^0 + ^^3(9^0 + ^) and p a _ - 

y y 

A'j = qiYoYsQj- In table 2 we gather our final continuum results for fi, in the MS scheme. For details 

about the analysis, results and its phenomenological implications we refer to [10] and [11]. 
5. D°-D° 

The bag parameters for the D° — D° oscillations can be determined following the same strategy 
as in [10] and outlined in the previous section, now with two of the OS valence quarks representing 
the up quark while the other two will be identified with the charm quark, i.e / ~ u and h ~ c. 
Physical values are obtained by interpolating data in /i» c " to the physical value pL c while chiral and 
continuum extrapolations are carried out simultaneously. The physical values for the charm quark 
mass have been previously derived and can be found in [13]. 

Exploratory studies show that using Gaussian smeared sources and choosing a time separation 
between meson sources smaller than T /2 (we set r sep = 16 at j8 = 3.8, T sep = 18 at jS = 3.9, 
?sep = 22 at j8 = 4.05) is crucial for quark masses around the physical charm and above. This 
improvement allows us to extract the ground state with more confidence and precision in a wider 
time interval. From figure la the benefit of Gaussian smearing compared to local source and sink 
is evident. For illustration in figure lb we display the quality plateau of Bj at the smallest lattice 
spacing, j3 = 4.05, and for the smallest value of the light quark. Preliminary results are collected 
in table 3. 



6. B°-t 

Since discretization errors on current lattices are expected to be large at the physical value of b- 
quark mass, our strategy for the computation of the bag-parameters in the B sector is based on the 



2 ^ sources are local for K-meson and smeared for D- and B- mesons 
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MS(2 GeV) 


Si 


B 2 


S 3 


B 4 


B 5 


0.78(04) 


0.71(04) 


1.10(12) 


0.94(06) 


1.17(15) 



Table 3: Continuum limit results for parameters of the D —D system renormalized in the MS 
scheme of [12] at 2GeV using Ml -type RCs defined in [1 1] 



ratio method approach proposed in [14] by introducing suitable ratios with exactly known static 
limit and interpolating them to the b-quark mass. 

A similar strategy to the one employed in [15] for the computation of the Bj an B s decay 
constants is followed here for the B B d and Bb s parameters. According to HQET, the ratio of renor- 
malized bag parameters evaluated in QCD is expected to approach unity as 1 / A/i — > 0. The leading 
deviation are predicted to be of order l/log^/Agco) and are expected to be tiny in the A/i range 
of our data. We consider the following chiral and continuum extrapolated scalar ratios, where the 
corrections to the power scaling in 1 / A/i (which are found to be at most ~ one standard deviation, 
see Eq.(6.5)) are estimated in perturbation theory by matching HQET to QCD through the C-factors 

<°d (rf ) = lim Aiefl ^ Aa/rf lim^o<0rf (A W ;A^,A/,a) , 



Ms (A/? ) = lim/i^a,,, lim a ^o fi£ (A {n) ; /Ui , fls , a) 



(6.1) 

C (A^ A*,a0 B Bd/s (f^ n) ; kea, A/„fl) 

C [fll VA;A*,AtJ Bb<i/s (A* 7*;A«»,A//*>aJ 
The C-factors ratio contains the information on the l/log(A/i) corrections. At tree-level (TL) 

C ( A/j' !) I A*> M ) = 1» while at leading-log (LL) it is given by 



/ ^ (n ) ^ * 

c Aa ;m ,m 



a(A*)/«rf ) ) «(aO/«(m) a(A*)/«« ; ) 



>) 



-(7^)/2)8o 



7a /A) 



and at NLL (not included in this preliminary analysis) the HQET mixing of 0\ with 02,3 should 
be considered. In the formula above % and 1 characterize the anomalous dimension (AD) of the 
axial ^/-current and the "1 1" element of the AD-matrix of the hlhl operator in HQET. 

In order to have better control on the chiral extrapolation, we consider the double ratio Bbs/BbcI 

rL(»M , , \ (A w ;A^,A.v,a) 

Co \l*> W \fLna,fk,a) = L ), (n) ; ; r (6.3) 

which also tends to 1 in the continuum and chiral limit. The ratio method is also applied to the phe- 
nomenologically interesting quantity E, = (/g s / fsd)\/B Bs / B Bd by forming the ratios analogous 
to Eq.(6.3) at successive values of the heavy quark mass. 

The quantities 00q, cof, ^ and Q have a smooth chiral and continuum extrapolation, showing 
no significant dependence on \i\ and small cutoff effects. The results turn out to be well described 
by a linear dependence in \ii and a 2 . For instance in figure 2a we show the chiral and continuum 
extrapolation of at the fourth of the considered masses. 

Finally, we study the dependence of the ratios (Od, (O s , Cm and £e on the inverse of the heavy 
quark mass as shown in figures 3a, 3b and 3c for a s , and Q. For co d and co s we perform a linear 
{ m d/s = 1 +t>(X)/fih) or a quadratic (aw, = 1 +b(X)/fif l + c(X)/fif J ) interpolation to the b-quark 
mass and the difference between them will be eventually assigned as a systematic error. In contrast, 
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£o and Q show very weak dependence on the heavy quark mass, thus in these cases we perform 
either a quadratic or linear interpolation or we fix the ratio equal to its asymptotic heavy-quark 
mass limit £ ffl = 1 and £e = 1. 

From the structure of Eq.(6.2) we derive the iterative formulae 



,V.(2)\ /-(3)\ ( \(K+1) 



C«(A? , )CefAf , )-CefAf +1) 



c(^)^(Af +1) 
c(Ai 1) )^(tf ) 

Bbs/Bbj (Af +1 

(a! 

fBslfBdy/Bas/BBd +1) 



/W /firf y/BBs/Bsd ( Aft 



relating the quantity at the heavy quark mass A/T with its triggering value at A{ and the fit 
values of ft) f /, ft) s , and £*, where the value of A/f together with the A factor and the number 
of steps to arrive to the b quark mass point have been previously tuned in such a way that after a 
finite number of steps (K) the heavy-light meson mass assumes the experimental value of Mb [17]. 
The value of the quantities in Eq.(6.4) at the triggering mass point can be accurately measured in 
the chiral and continuum limit since lies in the well accessible charm quark mass region. As 
example, in figure 2b we display the continuum and chiral extrapolation of the ratio Bss/Bgd at trig- 
gering mass, which shows small cutoff effects. Finally, Eq.(6.4) leads to the following preliminary 
estimates: 

5^ = 0.87(05), B Bs = 0.90(05) (6.5) 

where the error is the sum in quadrature of statistical and systematic uncertainties coming from the 
variation of the result if we use a linear or quadratic fit in the heavy quark mass and the truncation 
of the C-factors estimated as the difference between TL and LL matching. C-factors cancel out in 
the ratios Bb s /B b ^ and while the implementation of a constant, a linear or a quadratic fit in A/i 
turns out to provide essentially identical results. The preliminary estimates are 

B Bs /B Bd = 1.03(2), § = 1.21(06) (6.6) 

where in the second one we added in quadrature a systematic error due to the fit ansatz at the 
triggering mass point. 
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(a) (b) 
Figure 2: Double ratio for the fourth analysed mass (a) Chiral and continuum extrapolation of 

the ratio Bb s /Bb c i computed at the triggering point (b). Vertical lines represent the position of the 

physical point 
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